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Abstract 

Increasing interest is being dedicated in the last few years to the issues of exact 
computations and asymptotics of spin networks. The large-entries regimes (semiclas- 
sical limits) occur in many areas of physics and chemistry, and in particular in dis- 
cretization algorithms of applied quantum mechanics. Here we extend recent work 
on the basic building block of spin networks, namely the Wigner 6j symbol or Racah 
coefficient, enlightening the insight gained by exploiting its self-dual properties and 
studying it as a function of two (discrete) variables. This arises from its original defi- 
nition as an (orthogonal) angular momentum recoupling matrix. Progress also derives 
from recognizing its role in the foundation of the modern theory of classical orthog- 
onal polynomials, as extended to include discrete variables. Features of the imaging 
of various regimes of these orthonormal matrices are made explicit by computational 
advances -based on traditional and new recurrence relations- which allow an interpre- 
tation of the observed behaviors in terms of an underlying Hamiltonian formulation as 
well. This paper provides a contribution to the understanding of the transition between 
two extreme modes of the 6j, corresponding to the nearly classical and the fully quan- 
tum regimes, by studying the boundary lines (caustics) in the plane of the two matrix 
labels. This analysis marks the evolution of the turning points of relevance for the semi- 
classical regimes and puts on stage an unexpected key role of the Regge symmetries of 
the 6j. 
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1 Introduction 



This is an account of progress on understanding elementary spin networks in view of 
their ubiquitous occurrence in applications far beyond the theory of angular momentum 
in quantum mechanics where they were introduced by Wigner, Racah and others. They 
have been (and are) precious tools for computational nuclear, atomic, molecular and 
chemical physics: best known examples are vector coupling and recoupling coefficients 
and 3nj symbols CD E] EJ ffl EJ □ E E M HU El El [HE]- 

The diagrammatic tools for "spin networks" were developed by the Yutsis school 
and by others lfT6l[T7l . and were given this collective name in connection with applica- 
tions to discretized models for quantum gravity after Penrose [18 ], Ponzano and Regge 
fl9l , and many others (see e.g. I20,12"T1). 

The basic building blocks of (all) spin networks are the Wigner 6j symbols or 
Racah coefficients, which we are study here by exploiting their self dual properties 
and consequently looking at them as functions of two variables. This approach is most 
natural in view of their origin as matrix elements describing recoupling between al- 
ternative angular momentum binary coupling schemes, or between alternative hyper- 
spherical harmonics, or between alternative atomic and molecular orbitals, etc., and 
utilizes progress in understanding their role in the foundation of the modem theory of 
classical orthogonal polynomials, as extended to include discrete variables |22|. Fea- 
tures of the imaging of the orthonormal matrices is made possible by computational 
advances, that permit to elaborate accurate illustrations and comparisons, using exact 
computations based on traditional and new recurrence relations, which allow in turn an 
interpretation of the observed behaviors in terms of an underlying Hamiltonian formu- 
lation. An unexpected key role of the (mysterious) Regge symmetries 11231 of the 6j 
is briefly discussed. Suitable use of the results for discretization algorithms of applied 
quantum mechanics is stressed, with particular attention to problems arising in atomic 
and molecular physics. 

For background information and notation, we will refer to our previous papers 
ll7ll9l [TTIl . regarding the one dimensional view, and also to ll24l l25ll . where ample at- 
tention is devoted to the two dimensional perspective and where the related 4j model 
is elaborated in great detail. 

2 Background information 

Semiclassical and asymptotic views are introduced to describe the dependence on pa- 
rameters. They originated from the early association due to Racah and Wigner to geo- 
metrical features, respectively a dihedral angle and the volume of an associated tetra- 
hedron (see Fig. 4 in ifTTl ). which is the starting point of the seminal paper by Ponzano 
and Regge [ 19 1. Their results provided an impressive insight into the functional depen- 
dence of angular momentum functions showing a quantum mechanical picture in terms 
of formulas which describe classical and non-classical discrete wavelike regimes, as 
well as the transition between them. 

The symmetries with respect to exchange of role of the matrix entries of the Qj 
(transpositions), easily understood also from the associated tetrahedron picture, are 
well characterized. Similarly, we find that the mirror symmetries can be exploited 
when the symbols are employed in applications beyond angular momentum theory, 
implying meaning of the entries as quantities that can have a sign. 
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From the detailed study of combinatorial properties (requiring only triangular re- 
lationships associated to angular momentum coupling, and the closure quadrilateral 
property associated to coupling of four angular momenta), it emerged very recently 
1 24 1 how to uniquely assign labels to the elements of the grids for the matrix for which 
the orthonormalized Racah-Wigner coefficients are the elements. Figures of zero vol- 
ume and ridges are among properties that can be monitored on such square "screen": 
we sketch relationships and concept of Regge twins, "canonical" form and simplifica- 
tion of the following due to a suitable chosen ordering. From a computational view- 
point, explicit formulas are available as sums over a single variable. However, resorting 
to recursion formulas appears most convenient for exact calculations. Also, we will 
exploit them for semi-classical analysis, both to understand the high j-limit and, in 
reverse, to interpret the symbols them as discrete wave-functions obeying Schrodinger 
type of difference (rather then differential) equations. We have derived and computa- 
tionally implemented a two-variable recurrence that permits construction of the whole 
orthonormal matrix. The derivation follows our paper in 1 10 1 and is of interest also for 
other 3nj symbols in general. Separation of the two-variable recurrence relation leads 
to the basic three-term recurrence as depending on a single variable. It can be shown 
that it can be re-derived also from the Biedenharn-Elliot relationship in a form that 
shows the connection to a Schrodinger type of equation in the Hamiltonian formula- 
tion (for alternative Lagrangian formulations, see [25]). The following analysis of the 
caustics (and of the ridge curves, see below) is intended to characterize the modes of 
the spin network, as well as the guiding principles of the asymptotic analysis. 



3 The screen: mirror, Piero and Regge symmetries 



Following e.g. the Schulten and Gordon approach |26|, in |27| it is shown that the 6j 
symbol becomes the eigenfunction of the Schrodinger-like equation in the variable q, 
a continuous generalization of ji 2 : 



d 



*(<z) = o 



(1) 



(2) 



where ^f(q) is related to 

ji h J12 

33 3 323 

and p 2 is related with the square of the volume V of the associated tetrahedron (Fig 4 
in ifTTI '). whose edges are considered continuous and given by J\ = j\ + 1/2, J 2 = 
j 2 + 1/2, J 12 = j 12 + 1/2, J 3 = j 3 + 1/2, J = j + 1/2, and J 23 = j 23 + 1/2. The 
Cayley-Menger determinant permits to calculate the square of the volume of a generic 
tetrahedron in terms of (squares of) its edge lengths according to 
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(3) 



The condition for the tetrahedron with fixed edge lengths to exist as a polyhedron 
in Euclidean 3-space amounts to require V 2 > 0, while the V 2 — and V 2 < 
cases were associated by Ponzano and Regge to "flat" and nonclassical tetrahedral 
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configurations respectively. Equivalent to ([3]) is the Gramian determinant, used in 11241 
and (25 1, which embodies a clearer relationship with a vectorial picture. 

Major insight is provided by plotting both Qjs and geometrical functions (volumes, 
products of face areas, etc.) of the associated tetrahedra in a 2-dimensional j'12 — j'23 
plane (the square "screen" of allowed ranges of jxi an d J23 to be used in all the pictures 
below). Actually both (|3]l and the Gramian are equivalent to the famous formula known 
to Euler but first found five centuries ago by the Renaissance mathematician, architect 
and painter Piero della Francesca. The formulas of course embody all the well known 
"classical" symmetries of geometrical tetrahedra, which show up in the 6j symbol as 
well on applying concerted interchanges of its entries. Non-obvious symmetries of 
particular relevance in what follows are listed below. 

(i) The mirror symmetry. The appearance of squares of tetrahedron edges entails 
that the invariance with respect to the exchange Jo — J implies formally 
3 ~ 3 ~ 1 w i m respect to the entries of the 6j symbol. Although this is 
physically irrelevant when the js are pseudo-vectors, such as physical spins or 
orbital angular momenta, it can be of interest for other (e.g. discrete algorithms) 
applications. Regarding the screen, it can be seen that actually by continuation of 
the abscissa x — J12 and ordinate y = J 23 to negative values, one can have repli- 
cas that can be glued by cutting out regions shaded in |9|. This allows mapping 
of the screen to the S 2 phase space found in [ 25 1 . 

(ii) Piero line. In general, an exchange of opposite edges of a tetrahedron (or of the 
two entries in a column in the 6j symbol) corresponds to different tetrahedra and 
different symbols. In Piero formula, there is a term due to this difference that 
vanishes when any pair of opposite edges are equal. In general, a line can be 
drawn on the screen when ranges of jx% and j'23 overlap and in the screen one 
may encounter what we call a Piero line: when two entries in a column are equal, 
this line is a diagonal corresponding to the exchange x «-» y. Then, images on 
the screen will be symmetric with respect to such a line, as shown below with 
explicit examples. 

(iii) Regge symmetries. The manifestation of these intriguing symmetries in the 
present context is of paramount interest: it might be elucidated through con- 
nection with the projective geometry of the elementary quantum of space, which 
is being reconsidered from the viewpoint of association to the polygonal inequal- 
ities (triangular and quadrilateral in the 6j case), which have to be enforced in 
any spin networks. We find insightful to exhibit the basic Regge symmetry ll23l 
in the following form 



ji h 312 1 = f (ii + h -33 +i)/2 (ji +h +h - j)/2 ju 

33 3 323 J 1 (-ji + h +33 +i)/2 (ji -h +33 +i)/2 J23 

31 + P 32 - P 312 
33 + P 3 - P 323 



(4) 



where p = (-ji + j 2 - j 3 - j)/2 = [(j 2 + j) - (ji + ia)]/2. Often the Regge 
symmetry is written in terms of the semi-perimeter s = (ji +j 2 + J3 +j) /2. Obviously, 

s = p + 3x+ 33 =~P + h + j- 

The key observation is made in ESI that the range of both J\ 2 and J23,namely the 
"size" of the screen, is given by 2min{ Jj, J 2 ,Jz, J, J\ + p, J 2 + p, J3 + p, J + p). 
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Therefore, for definiteness and no loss of generality, we are allowed in most cases to 
conveniently choose for the discussion whichever of the two 6j symbols twinned by 
Q contains such a minimum value. 



4 Features of the tetrahedron volume function 

Looking at the volume V as a function of x = J12 and y = J23, and after some alge- 
braic manipulations, we get the expressions for the xymax and yvmax that correspond 
to the maximum of the volume for a fixed value of x or y (respectively), namely 



A(J, J2) + J 2 y 2 — y 4 
2y 2 

A( J 2 , J) + J? x 2 - x 4 



■^Vmax o „,2 ' ^ ^ 



yVmax 2 j.2 ' ^ 



where 



A(a,b) = (Ji + a)(Ji-a)(J 3 + b)(J 3 -b), 

J? = Jf + Jl + Jl + J 2 - (7) 

We call "ridge" curves the plots of Eqs. |5]) and ([6]) on the x, y-screen. Each one marks 
configurations of the associated tetrahedron when two specific pairs of triangular faces 
are orthogonal. The corresponding values of the volume (V ma x,x an d V maXyy ) are 

v 2F(J 1: J,y)F(J 2 ,J 3 ,y) 
„ _ 2F(J 1 ,J 2 ,x)F(J,J 3 ,x) 

Vmax.y — ) K?) 

6X 

where 

^/(a + 6 + c)(a + & - c)(a - 6 + c)(-a + b + c") 
F(a, fe, c) = , (10) 

is the area of the triangle with sides a, b and c. 

Curves corresponding to V = (the caustic curves) obey the equations 



J V max 



12 V„. 



il 



Figure 



1(a) 



1 9 V 

2 2 v max,y „ n 
i/* = ^max± "> (11) 

shows these curves for the < _„ ^ "( 12 i case. The inside region 

t 55 60 j 23 J 

enclosed by the ovaloid is that of finite volume tetrahedra, the ovaloids themselves 
correspond to configurations of flattened tetrahedra, specifically convex planar quadri- 
laterals at the upper right corner,concave planar quadrilaterals at both the upper left and 
and the lower right corners, and crossed planar quadrilaterals at the lower left corner. 

For particular values of j%, j 3 and j linear configurations are allowed as well. 
These interesting degenerate cases are illustrated in the other panels of Fig. [T] In Fig. 



1(b) the case of j\ +j = j 2 +j 3 is considered. The other interesting cases are obtained 



for ii + 32 = h + J (Fig- 1(c) i and for j x + j 3 = j 2 + j (Fig. 1(d) i 
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Figure 1: Plots of the caustics and ridges given by Eqs.(5][TTi for four sets of 6j- 
symbols. Panel [1(a)] jt = 45, j 2 = 30, j 3 = 55 and j = 60. 15 < J 12 < 76 
and 25 < J 23 < 86. Panel [KE)] ji = 140, j 2 = 130, j 3 = 110 and j = 100. 
10 < J12 < 211 and 40 < J 23 < 241. Panel [T(c)l j 1 = 140, j 2 = 100, j 3 = 110 and 
j = 130. 40 < J12 < 241 and 10 < J 23 < 211. Panel j\ = 140, j 2 = 110, 
j 3 = 100 and j = 130. 30 < J 12 < 231 and 10 < J 23 < 211. 



Note that, due to the symmetry properties of the 6j symbols, the cases j\ + j = 
32 + j 3 an d ji + 3i — j 3 + j are intrinsically equivalent. Figure 1(d) shows the effect 
of degeneracy with respect to the far from obvious Regge symmetry, which manifests 
as specularity with respect to the Piero line, in this case the diagonal from the origin in 
the square screen. 



5 Symmetric and limiting cases 

When some or all the j's are equal, interesting features appear in the screen. Similarly 
when some are larger than others. 
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Figure 2: Plots of the caustics and ridges (Eqsf5]|TT[> for the 6j-symbols with j% = 100, 
h = 100, js = 150 and j = 210. 60 < J u < 201 and 110 < J 23 < 251. 



5.1 Symmetric cases 

For ji = j 2 plots of Eqs.([5pT|i like those of Fig. [2] are obtained (equivalent to the 
js = j one, in virtue of standard symmetries). 

For ji = j'3 plots of Eqs.( 5JTT i like that of Fig. [3]are obtained. For symmetry this 
case is equivalent to the ji = j one. 

Imposing bot h ji + j = j 2 + j 3 and j x + j 2 = js + j we have that j x = js and 
j 2 =j. See Figs. |4(a)1e|4(b) 

Imposing both j x + j = j 2 + js and j\ + j 3 = j 2 + j we have that j x = j 2 and 
js — j- In these conditions, plots of Eqs.( 5JTT i like those of Fig. |5]are obtained. 

This case is formally equivalent to the one where j\ = j and js = j 2 which is 
obtained imposing jx + j 2 = j 3 + j and j\ + j 3 = j 2 + j. 



5.2 Limiting cases 

Interestingly, the Fig. [5]permits us to discuss the caustics of the 3j symbols as limiting 
case of the corresponding Qj where three entries are larger than the other ones, namely: 
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100 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ■— »J 

100 120 140 160 180 200 220 240 

Figure 3: Plots of the caustics and ridges (Eqs|5fTT)i for the 6j-symbols with j\ = 100, 
j 2 = 150, j 3 = 100 and j = 210. 110 < J 12 < 251 and 110 < J 23 < 251. 
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Figure 4: Plots of the caustics and ridges (Eqs|5p"T) for two sets of 6j symbols. Panel 
|3(a)] ji = 200, j 2 = 100, j 3 = 200 and j = 100. 100 < J 12 < 301 and 100 < J 23 < 
301. Panel [4(b)] ji = 110, j 2 = 100, j 3 = 110 and j = 100. 10 < J 12 < 211 and 
10 < J 23 < 211. 



( 31 h h ) f Ji Ja js 1 n2) 

^ m x m 2 m 3 ^ j£5S„ \ h+R h + R l 3 + R J ' U ' 

where 

mi = Z 3 - l 2 , 

m 2 = Zi - Z 3 . (13) 

If we define 

mi~F + D, 
m 2 := F - D , 
R:= R + l 3 - D , (14) 

we have 

^ mi m 2 m 3 J R + F R-F R + D J v; 

For # oo it is ± F ~ 

The caustic of the 3j symbol is defined as 
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Figure 5: Plots of the caustics and ridges (Eqs|5fTT)i for the Qj symbols with j\ = 1000, 
j 2 = 1000, j 3 = 100 and j = 100. < J 12 < 201 and 900 < J 23 < 1001. According 
to the text, this figure also exemplifies features of ridges and caustics for a 3j symbol. 



10 



and 




(17) 

The case of caustics on the screen for four large entries, leading to reduced Wigner 
d matrices, is studied in 0. 

For ji = j2 = js = j plots of Eqs.pfTT) like those of Fig. [6] are obtained. This 
case is obtained when ji + j = j 2 + J3, ji + h = h + j and j x + j 3 = j 2 + j. 

Fig. [7] shows the caustics of a full-symmetric 6j when the symmetry is broken by 
varying j. 

6 Remarks and Conclusion 

Explicit computational formulas are available either as sums over a single variable and 
series. However recourse to recursion formulas appears most convenient for both ex- 
act calculations and -as we will emphasize- also for semiclassical analysis, in order 
to understand high j limits and in reverse to interpret them as discrete wavefunctions 
obeying Schrodinger type of difference (rather then differential) equations. In further 
work, we derive and computationally implement the two variable recurrence that per- 
mits construction of the whole orthonormal matrix The derivation follows our paper in 
1 10 1 and is of interest also for other 3nj symbols in general. 

The extensive images of the exactly calculated 6j's on the square screens illustrate 
how the caustic curves studied in this paper separate the classical and nonclassical 
regions, where they show wavelike and evanescent behaviour respectively. Limiting 
cases, and in particular those referring to 3j and Wigner's d matrix elements can be 
analogously depicted and discussed. Interesting also are the ridge lines, which separate 
the images in the screen tending to qualitatively different flattening of the quadrilateral, 
namely convex in the upper right region, concave in the upper left and lower right ones, 
and crossed in the lower left region. 

Interestingly, the Regge symmetries not only restrict the range of the discrete man- 
ifold, but also allow an assignment of the involved modes. 

As a continuation of our recent work, this paper has demonstrated that interesting 
insight into properties of the basic building blocks of spin networks, the Wigner 6j 
symbols or Racah coefficients, is gained by exploiting their self dual properties and 
studying them as a function of two variables, an approach most natural in view of 
their origin as matrix elements. Features of their imaging of the orthonormal matrices 
are fostered by computational advances, which is being developed on traditional and 
new recurrence relations, which also allow interpretation of the underlying Hamilto- 
nian mechanics: the borderline of the two limiting modes -quasi-classical and deeply 
quantum- has been the object of this paper. The characteristic features of this boundary 
line, described by the caustic curves which follow the turning points of the classical 
motion, have been elucidated. A key role was also revealed of the surprising Regge 
symmetries. 
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Figure 6: Plots of the caustics and ridges (Eqs|5pT[) for the full symmetric 6j-symbols 
with ji = 100, j 2 = 100, j 3 = 100 and j = 100. < J 12 , J 23 < 201. This figure 
clearly points out the interest of extending the screen by mirror symmetries. 
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Further work based on these tools involves detailed studies of some of the important 
aspects considered too briefly or not at all in this paper, such as the classical, mirror and 
Regge symmetries, the limiting cases to the simpler 3j coefficients and rotation matrix 
elements, the extensions to higher 3nj symbols and to the so called q deformations. 

Semiclassical and asymptotic analysis provide limiting relationships converging 
into the Askey scheme. Relationships also arose then with orthogonal polynomials, 
and indicate avenues to views to generalizations (continuous extensions, relationships 
with harmonics of rotation groups, and finally q-extensions, . . .). All these extensions 
may need some modifications when detailed properties are discussed, but in general 
studying 6j symbols continues to deserve further work. 

Use is finally pointed out for discretization algorithms of applied quantum mechan- 
ics, particular attention being devoted to problems encountered in atomic and molecular 
physics. 
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